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In this paper, we extract from concurrence its variable part, denoted A, and use A as a time- 
dependent measure of "distance" , either postive or negative, from the separability boundary. We use 
it to investigate entanglement dynamics of two isolated but initially entangled qubits, each coupled 
to its own environment. 
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I. INTRODUCTION 



Entanglement is arguably the most intriguing feature 
of quantum mechanics [l| . In a ideal world, the entangle- 
ment needed for quantum information processing would 
be stable and uncorrupted. But in reality, the interaction 
of a quantum system with its surroundings is unavoid- 
able because no real physical system can be isolated com- 
pletely from external noise. The effect of noise on pairs 
of quantum systems is particularly interesting when the 
systems are relatively well isolated from noise and un- 
dergo relaxation individually rather slowly, but may be 
vulnerable to qualitatively different relaxation channels 
in regard to mutual coheence, i.e., to entanglement. Five 
years ago such a difference was pointed out by the group 
of Knight 0. The present paper can be considered an 
indirect extension of the results reported there. 

More recently new features have been discovered for 
entanglement decoherence as a time-dependent process. 
It has begun to be the subject of wide interest H, 0, H, 
I, 0, S, H M M M M m EE QE E3- Contrary to 
intuition based on experience with single qubit decoher- 
ence, we have shown [5] that entanglement may decay to 
zero in a finite time (and remain zero for at least a finite 
time). 

This finite-time decay is called entanglement sudden 
death (ESD) . It is a feature of nonlocal quantum coher- 
ence, one that introduces a new bipartite element into 
quantum open system behavior [IE, HE] • ESD adds an el- 
ement of complexity to the commonly encountered state- 
space diagrams showing evolution along various routes 
that, in the normal cases of relaxation lead to separabil- 
ity. An example is shown in Fig. [TJ where trajectories 
all begin as entangled states (in the non-separable space) 
and evolve toward separability. 

Here we identify the boundary between the S and NS 
regions as a border that establishes the origin of a gen- 
eralized "coordinate" , which we define. In this context 
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FIG. 1: Quantum trajectories in a non-metric "evolution space" 
of density matrices of two qubits. The space of the two qubits can 
be divided into the region of separable states denoted by S, and 
entangled states denoted by NS. Within the separable states S we 
identify an important subset (ESD) that includes all the states that 
reach separability abruptly in a finite time rather than asymptoti- 
cally. The diagram shows the three initially entangled states pi , p2 
and P3 evolving into three separable states. 



we examine the evolution of different initial NS states 
under different decoherence mechanisms as their trajec- 
tories pass the boundary, i.e., as their coordinate distance 
from the boundary passes from negative to positive val- 
ues. We show situations in which this border-crossing 
behavior is (a) monotonic or (b) periodically repeating, 
or (c) what we think of as the most familiar type of de- 
coherence evolution, reaching the border but not crossng 
it. 



II. THE A DISTANCE FOR PURE AND MIXED 
STATES 



For a two-qubit state described by density matrix p, 
Wootters' concurrence C{p) is the most-used measure 
of entanglement (non-separability) because of its rela- 
tive ease of direct calculation [19(, and its conveniently 
normalized range from C = for a separable state to 
C = 1 for a maximally entangled state. It should be 
emphasized that the boundary of separability of a state 
is independent of the entanglement measure, and in this 
respect concurrence is exactly in accord with entropy of 
formation, negativity and tangle, while many other fa- 
miliar coherence concepts such as purity or non-zero di- 
agonal matrix elements have no reliable relation to en- 
tanglement. 
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For two qubits, the concurrence may be calculated ex- 
plicitly from the density matrix p for qubits A and B: 



= max (0, - yfi* - - yfc) , (1) 

where the quantities A^ are the eigenvalues in decreasing 
order of the matrix 



C = P«®<)P*«®<), (2) 
where p* denotes the complex conjugation of p in the 



j ■ / ' i ' /'i ' 1 / > i ' 



Pauli matrix expressed in the same basis as: 



A,B 







(3) 



The relation among these eigenvalues is of course the 
essence of the concurrence, so for convenience we define 
the quantity A: 

A = ^/\[ - y% - y% - V^4- ( 4 ) 

It is obvious that A can serve as a one-dimensional quan- 
tum measure of the "distance" of a state away from the 
boundary and whether it is on the positive (entangled) 
or negative (disentangled) side of the solution space of 

Fig. m 

A focus on A is prompted because it contains differ- 
ent information than the concurrence C (p) from which it 
originates. To see this and to begin our discussion, we 
evaluate A for an arbitrary pure state that is separable 
(not entangled): p AB = \^f AB ){9 AB \, where 



b\-) 



d\-) 



It is easy to check that both C(p) = and the distance 
A = for this case. 

From these comments we show that distinctions be- 
tween C and A make A more useful. For example, we 
see that, given any density matrix p, if A(p) is strictly 
negative, A < rather than A < 0, it implies that p 
must be both mixed and separable. Note that this is 
more than the value of C(p) provides. It cannot tell if 
the state p is mixed or not. Another unusual property 
of A during decoherence evolution is that its final value 
will generally depend on the system's initial state, even 
under fully disentangling evolution, whereas a fully dis- 
entangled concurrence must have the value no matter 
what the initial state was. From the final value of A one 
can determine if ESD has occurred, but not from the final 
value of concurrence. 

Clearly it is important to be aware if a boundary- 
crossing transition occurs for a bipartite A — B unit in 
an active quantum network; and the existence of entan- 
glement sudden death (ESD) signals that the transition 
may occur abruptly and in that sense very early, as ex- 
perimental work is beginning to demonstrate [20(. Here 
we describe the way ESD relates to changes in the A dis- 
tance for several state trajectories. In all cases we will 
suppose that A and B are initially entangled and experi- 
ence relaxation that drives the pair to a separable state. 



III. PHASE-RELAXING EVOLUTION OF A 

First we consider the time evolution of the A measure 
under broadband dephasing. The following familiar spin 
Hamiltonian 0, i, El provides a good illustration. We 
take h — 1 and write: 



H(t) = --n(b A (t)af + b B (t)a l z 



(5) 



where p stands for a common gyromagnetic ratio and the 
magnetic fields b A (t) and bsit) are assumed to be weak, 
independent, statistically stationary and noisy. We as- 
sume their bandwidths are sufficiently smooth and broad- 
band in nature to exclude the existence of decoherence- 
free subspace effects. They have the average values: 



< bi(t) > = 0, 



< bi(t)bj(t') > 



^5 lj S(t-t'), i,j = A, B, 



(0) 
(7) 



where the Ti are the phase damping rates of qubits A and 
B. For this interaction, the time dependence of density 
matrix p AB (t) of the two qubits has been found [3| and 
is given by: 



p AB (t) = 

( Pll 7-BP12 

7-BP21 P22 

7AP31 JA1BP32 

\jA-fBP41 7AP42 



7-AP13 7^7B/°14\ 
1A1BP1Z JAP24 ] 



P33 
7BP43 



lBP3i I 
P44 / 



(8) 



where the pij elements here are initial values and the 7's 
contain the time dependences: 



1A 



-t/2Ti 



IB 



-t/2T 2 E 



0) 



For convenience we will take the relaxation times to be 
equal: T% = T 2 B = 1/T. 

For our purpose it is enough to examine the evolution 
of A for an initial state in the more specialized class of 
bipartite "X" density matrices 0: 



P AB (0) 



a w\ 

b z 

z* c 

\w* dj 



(10) 



where a + b + c + d= 1. The class of mixed X states 
arises naturally in a wide variety of physical situations, 
which include pure Bell states as well as the well-known 
Werner mixed states [21|. Here we will take w = 0. 

From the general solution ([5]), one can easily show for 
the initial state ([Toll that one finds for w = 



(11) 



A(t) = 2|z(i)| - 2\ 



The diagonal elements are independent of i, but since z(t) 
goes to zero as exp(-Tt) as t increases, A must become 
strictly negative if ad ^ 0. This condition mandates that 
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FIG. 2: Two graphs of state trajectories are shown. On the left 
two A curves leading to a final separable state, and to the same final 
value of concurrence, do not reach the same A value. The "familiar" 
asymptotically smooth relaxation of the solid curve, reaching full 
separability only in the limit t — > oo, is obtained with a = 0, b = 
c = d = z = 1/3. The dashed curve, obtained for a = d = 
1/12, b = c = z = 5/12, exhibits ESD by reaching separability 
at a finite time when A becomes negative. In the first case A 
never becomes negative. The non-metric evolution space picture 
of the same processes is shown on the right, emphasizing that one 
trajectory crosses the boundary and the other only asmptotically 
touches it. 

ESD occurs. The contrasting results that are possible 
are shown in Fig. [2] One can also check that if both 
ad and be are nonzero, then ESD also always occurs for 
the initial mixed density matrix (|10[) with both w and z 
non-zero. 



types of Bell states: 

\$ab) = cosa|+,+) +sina|-,->, (13) 

and 

\^ab) =cosor|+,-) +sma|-,+). (14) 

A surrogate for "reservoir" action here is introduced by 
tracing the radiation modes. This acts as a decoherence 
process for the A— B pair because it discards completely 
any information stored in knowledge of the radiation, and 
as in the preceding example leaves A — B separability as 
the issue at hand. An earlier examination in which single- 
mode relaxation was studied [2J included an additional 
physical element by permitting thermal statistics for the 
mode. In our simpler case, for the initial state (| 1 3D the 
time-dependent expression for A turns out to be [lq : 

A(i) = cos 2 (Gt/2)^2| cos a sin a\ - 2 sin 2 (Gi/2)| cos of) . 

(15) 

Both of these are shown in Fig. [3] 

In contrast, for the initial state (fl4)) . the resulting A is 
given by 

A(t) = |sin2a|cos 2 (Gi/2). (16) 

In the second JC case A > for all possible a and 
t, and A is not zero in any finite time interval (except 
for the trivial case where the initial state is separable). 
We immediately conclude that ESD never happens for 
this type of initial state. However, in the first case the 
function A(t) can take negative values. For example, for 
any a such that | cosa| > | sina|, as t increases from t — 
the function A(t) will take negative values at finite times 



IV. EVOLUTION OF A IN CAVITY QED 

In our second example each of a pair of two-level atoms 
is confined in a high-Q cavity and exposed to a single 
mode of radiation resonant with the atomic transition. 
Each separate atom-mode interaction produces a Jaynes- 
Cummings evolution governed by separate Hamiltonians 
[H]. With h= 1, we have 

#tot = + ( ^{a)a A + a J la)+uja)a 

+ ^ + ^a B +afb)+^b (12) 

where wo and u> are the frequencies of atoms and cavities, 
and G is the atom- mode coupling constant, the vacuum 
Rabi frequency. This model has four interacting parties 
that we can label in an obvious notation A, B, a, b, and 
it permits well-known exact solutions. 

Here our point is easily illustrated by taking exact res- 
onance and vacuum initial states for the modes and by 
examining just two sets of entangled initial states for the 
two atoms A and B, namely superpositions of the two 




FIG. 3: Two representations of the same decoherence processes 
are shown. On the left, two curves of A vs. t are shown, the first 
for initial state ^ab and the dashed curve for state <t>AS > in both 
cases for a = 7r/6. On the right these periodic processes are seen as 
periodic trajectories in evolution space, with the interesting feature 
that one of them repeatedly touches but never crosses the boundary 
to separability. 
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t, and also remain negative for finite lengths of time. 
These contrasting behaviors are shown in Fig. We 
see immediately striking differences, compared with Fig. 
[31 arising from well-known periodic behavior induced by 
the pure JC interaction. 

V. CONCLUDING COMMENTS 

In summary, we undertook to examine the disappear- 
ance of entanglment in a finite time (ESD) in relation to 
considerations suggested by the generic sketch in Fig. [Tj 
of trajectories of states in "evolution space" . We consid- 
ered bipartite examples that are explicit, solvable exactly, 
and simple rather than complicated, to be clear that con- 
sequences are not arising from complications but rather 
from intrinsic properties. We treated only decoherence 
processes that carry initially entangled states into sepa- 
rable states. 

To carry out our discussion we used obvious proper- 
ties of the separability "distance" A. We followed two 
exactly solvable models in which quantum entanglement 
decreases due to the interaction of qubits A and B with 
their environments, while they are not interacting with 
each other, and their environments are also not mutually 
interacting. In the second example, the "environment" 
consists of only one cavity mode for each atom, but one 
can say it acquires the status of environment when it 
is traced out, i.e., when it is used as an information loss 
channel. We have pointed out that the A "distance" con- 
tains information that concurrence does not. In particu- 



lar, if A(p) becomes negative, it implies that A describes 
a separable mixed state. 

A few further comments are in order: (1) The degree 
of entanglement in this paper is measured by the gen- 
eralized concurrence A which is valid for both pure and 
mixed states. For pure separable state, one always has 
A = 0. This means that the negative A's only happen for 
a mixed separable state. (2) As usual, we emphasize that 
entanglement sudden death is independent of the choice 
of measures of entanglement. (3) We note that while we 
have found interesting and unexpected features of A for 
two solvable models, it would be very interesting to ex- 
tend these results to more general quantum dynamical 
systems. 

A natural question remains for further study: Why 
does bipartite ESD happen for some initial states, but 
not others, even if both ultimately become separable? 
The question may be addressed on several levels and we 
believe that a satisfactory answer still does not exist. 



Acknowledgments 

With this small note we add our congratulations at 
what is probably only the midpoint of the research career 
of Peter Knight. We are among the many admirers of his 
contributions to quantum information and to quantum 
optics more generally over many years. Financial support 
has been received from US Army Research Office grant 
ARO 48422-PH. 



[1] A. Einstein, B. Podolsky and N. Rosen, Phys. Rev. 47, 
777 (1935). 

[2] M.S. Kim, J. Lee, D. Ahn and PL. Knight, Phys. Rev. 

A 65, 040101 (2002). 
[3] T. Yu and J.H. Eberly, Phys. Rev. B 66, 193306 (2002). 
[4] T. Yu and J.H. Eberly, Phys. Rev. B 68, 165322 (2003). 
[5] T. Yu and J. H. Eberly, Phys. Rev. Lett. 93, 140404 

(2004). 

[6] L. Diosi, in Irreversible Quantum Dynamics, edited by F. 
Benatti and R. Floreanini (Springer, New York, 2003), 
pp. 157-163. 

[7] P.J. Dodd and J.J. Halliwell, Phys. Rev. A 69, 052105 
(2004). 

[8] L. Jakobczyk and A. Jamroz, Phys. Lett. A 333, 35 
(2004). 

[9] T. Yu and J.H. Eberly, Quant. Inf. Comp. (in press); 

arXiv: |quant-ph/0503089[ 
[10] D. Tolkunov, V. Privman, and P. K. Aravind, Phys. 

Rev. A 71, 060308 (2005). 
[11] M. F. Santos, P. Milman, L. Davidovich and N. Zagury, 



Phys. Rev. A 73, 022313 (2006). 
[12] S. Shresta, C. Anastopoulos, A. Dragulescu, and B. L. 

Hu, Phys. Rev. A 71, 022109 (2005). 
[13] F. Mintert, A. R. R. Carvalho, M. Kus, A. Buchleitner, 

Phys. Rep. 415, 4, 207(2005). 
[14] M. Ban, J. Phys. A 39, 1927 (2006). 
[15] T. Yu and J. H. Eberly, Opt. Commun. 264, 393 (2006). 
[16] M. Yonac, T. Yu and J. H. Eberly, J. Phys. B 39, S621 

(2006). 

[17] Z. FicekandR. Tanas, Phys. Rev. A 74, 024304 (2006). 
[18] T. Yu and J. H. Eberly, Phys. Rev. Lett. 97, 140403 
(2006). 

[19] W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 
[20] M. P. Almeida, F. de Melo, M. Hor-Meyll, A. Salles, 

S. P. Walborn, P. H. Souto Ribeiro and L. Davidovich, 
|quant-ph/0701184l 
[21] R. Werner, Phys. Rev. A 40, 4277 (1989). 
[22] E.T. Jaynes and F.W. Cummings, Proc. IEEE, 51, 89 

(1963). 



